Exact approaches to charged particle motion in a time- dependent flux-driven ring 
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We consider a charged particle which is driven by a time-dependent flux threading a circular ring 
system. Various approaches including classical treatment, Fourier expansion method, time-evolution 
method, and Lewis- Riesenfeld method are used and compared to solve the time-dependent problem. 
By properly managing the boundary condition of the system, a time-dependent wave function of 
the charged particle can be obtained by using a non-Hermitian time-dependent invariant, which 
is a specific linear combination of initial angular-momentum and azimuthal-angle operators. The 
eigenfunction of the linear invariant can be realized as a Gaussian-type wave packet with a peak 
moving along the classical angular trajectory, while the distribution of the wave packet is determined 
by the ratio of the coefficient of the initial angle to that of the initial canonical angular momentum. 
In this topologically nontrivial system, we find that although the classical trajectory and angular 
momentum can determine the motion of the wave packet; however, the peak position is no longer 
an expectation value of the angle operator. Therefore, in such a system, the Ehrenfest theorem is 
not directly applicable. 

PACS numbers: 03.65.Ge, 03.65.Fd, 73.21.-b 
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I. INTRODUCTION 

Charged particle driven by a time-dependent pertur- 
bation in a quantum system is a nontrivial fundamental 
issue pi 0, 0, 0, IE IE 0- O ne can access the charged 
particle wave function by placing it in a quantum ring 
threaded by a time-dependent magnetic flux. The vector 
potential A(i) associated with the time-dependent flux 
times the charge q leading to a phase shift propor- 
tional to the number of flux quanta penetratingthe ring , 
this is known as Aharonov-Bohm (AB) effect 0, |E llOf . 
In adiabatic cyclic evolution, Berry first discovered 
that there exists a geometric phase. Later on, Aharonov 
and Anandan (AA) removed the adiabatic restriction to 
explore the geometric phase for any cyclic evolution • 
Time-dependent fields are also used to deal with field 
driven Zener tunne ling , in which nonadiabaticity plays a 
crucial role 0QG]|. 

In mesoscopic systems, a number of manifestations of 
the AB effect have been predicted and verified 0, 0, 
111. ITg| . On the other hand, Stern demonstrated that 
the Berry phase affects the particle motion in the ring 
similar to the AB effect, and a time-dependent Berry 
phase induces a motive force |20|. Recently, AA phases 
have been found to play the role of classical canonical 
actions and are conserved in the adiabatic evolution of 
noneigenstates [2lj . 

In the present work, we consider a noninteracting spin- 
less charged particle moving cyclically in a quantum ring 
in the presence of a time-dependent vector potential. 
Such a particle motion can be described by the time- 
dependent Schrodinger equation 



time-dependent vector potential A(t), given by 
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P - qA(t) 
L - qRA(t) 



(2) 



Here P = egPg is the canonical momentum operator with 
eg being the unit vector along the azimuthal angle 0; 
L = L z = (r x P) z is the canonical angular momentum 
operator in the z direction; / = mR 2 is the moment of 
inertia of the particle; A(t) = A(t)eg is the vector poten- 
tial; and R is the radius of the circular ring. This time- 
dependent dynamical problem can be solved by taking 
into account the Fourier expansion, time evolution oper- 
ator, and Lewis- Riesenfeld (LR) method p2l |2^|. 



II. A CLASSICAL TREATMENT 

We first analyze the time-dependent problem in a clas- 
sical manner. The time- varying magnetic flux induces 
an electric field E = Eeg such that E — —dA/dt. The 
charged particle thus obtain a kinetic momentum incre- 
ment during the time interval from to t, namely 

A Pc = A(mv) = m[v(t) - «(0)] = -q [A(t) - A(0)} , (3) 

where p c — mv is the kinetic momentum. It should be 
noted that both p c and qA are not conservative quanti- 
ties, while from Eq. (0 we see that the canonical mo- 
mentum P c is a constant of motion: 

P c (t) = mv(t) + qA{t) = mv(0) + qA{0) = P c (0). (4) 

Comparing the two identities in Eq. (J2J, we see that 
the result of Eq. © is equivalent to 



where the Hamiltonian H(t) is induced by an external 



Al c =I[w(t) -w(0)] 
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[*(*)- *(())], (5) 
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where l c = (r x p c ) 2 = Iio indicates the kinetic angular 
momentum, lu is the angular velocity, $ is the magnetic 
flux threading the ring, and the fact 

$(*) = 2nRA(t) (6) 

has been used. Also, from Eq.QJ we have 

Lc{t) = IS) + J^W = lc(0) + ±*(0) = L B (0). (7) 

Thus the canonical angular momentum L c — (r x P c ) z is 
also a constant of motion. 

Now we define the writhing number as 



n$ (t) 



$0 ' 



(8) 



where $o = h/q is & flux quantum. We also define 

L c = n h, l c (t) = n c (t)h, (9) 

then we have 

n c (t) = no - n$(t). (fO) 

All of these n's are reaZ numbers. 

Now the angular position of the driven particle is given 
by 



.(t) = O + 



n c (r)h 



= 8q+ uj t - 



7i$(r)7i 



dr. (11) 



Here 9q indicates the initial azimuthal angle; and u>q = 
w(0) = noh/I stands for the initial angular velocity. 
Below we denote the initial kinetic angular momentum 
lo = Z c (0) for simplicity. 



III. A FOURIER EXPANSION METHOD 

The simplest method for solving the time-dependent 
flux-driven problem is the Fourier expansion method. 
The first thing about the system we discuss is that the 
wave function satisfies the periodic boundary condition: 



ip(B,t) =ip(e + 2-K,t). 



(12) 



The most general form of ip for the present problem is 
thus written as 



Substituting Eq. I|13fl into Eq. (JTJ, we can find 

oo 

53 ihc n Ut)e me 



53 °nfn{t) 



(nh~qRA(t))\ mg 
21 



(14) 



n— — oo 



Solving Eq. (|14|) . after some procedures we obtain 



f n (t) = exp 



2Ih 



ih- qRA(t)} 2 dt\ , (15) 



and thus 



ih /"* 

ip(0,t) = 53 c « ex P "i — 777 / [ n — n *( r )] dr + inO 



n— — oo 



21 



(16) 



As an simple example, let us choose 

c n = TV exp [-(T 2 (n - n Q ) 2 - i6 (n - n Q )] , (17) 

where N indicates an appropriate normalization con- 
stant; cr, #o and no are real numbers, and n is an integer. 
Substituting Eq.JUJ into Ea.lfTB^. we get 



ip(9,t) = Nexp 



TL— — OO 



h Jo 21 



it 



dr + ino9 c (t) 



(18) 



53 exp \ -a 2 ( 1 + - ) (n - n ) 2 + in [6 - 9 c (t)} 



where T = 21a 2 jh. Applying the Poisson summation 
formula 



oo oo , „ c 

E /(«) = E (/_ 



f(x)e l27rnx dx (19) 



on the function 
f(x) = exp 



if 



-<7 2 ( 1 + - ) (x - n Q ) 2 + i(9- 6 c {t)) x 



(20) 



we can obtain an alternative expression 



i>(6,t) = N 

OO 

x 53 ex p 

n— — oo 



a 2 (l + f) 



exp 



h J 21 



dr 



(21) 



(0 - c (t) + 2717T)" 

4a 2 (1 + |) 



+ ino(9 + 2mr) 



^(6,t) = 53 c„/„(f)e 



i nO 



(13) 



n— — oo 



where the c„'s are appropriate coefficients to be deter- 
mined by the initial and the boundary conditions. 



IV. A TIME EVOLUTION METHOD 

In this section, we shall present how to get the general 
solution shown in the previous section in terms of the 
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time evolution operator U(t). The state \ijj(t)) is con- 
nected with the initial state |?/>(0)) through 

\m) = u(t)\m)- (22) 

and the wave function ip(9,t) is given by 

V(0,i) = (0|#(tMO)>, (23) 

where \6) is the 0-eigenket in the Schrodinger picture, 
which will be explained later. 

To begin with, we introduce the canonical commutator 

[0(O),L(O)]=ift. (24) 
From this identity, we have 

[0(f), L(t)} - UHt)[9(0), L(0)]U(t) = i%. (25) 
Utilizing Ea. (|25|l we can derive 
dL{t) [L(t),H(t)] 



dt 



if i 



= 0. 



(26) 



we thus obtain the identity L(t) = L(0). Following similar 
procedure it is easy to obtain 



d9(t) [0(t),H(t)] L(0) - n$(t)n 



dt 

which gives us 



it) 



1 ./n J- 



(27) 



(28) 



We can see that the canonical angular momentum is a 
constant of motion. This is consistent with the classical 
results discussed in Sec. ITT1 

From the above results we have 



[H(t),H(t')}=0 



(29) 



for any two times t and t' . Hence the time evolution 
operator is simply given by 



U{t) 



exp 



exp 



ft Jo 



ih 



27. / V k 



H(r)dT 

(Up) 



n$(r) rfr 



(30) 



To proceed further, we define \n) as the eigenket of 
L(0) obeying 



L(0)\n) = ri%). 



(31) 



In addition, we also assume that \9) is an eigenket of 



J9(0) 



obeying 



oi<»(0) 



\6)=e*\d) 



(32) 



The orthogonal conditions of the two eigenkets are given 
by 



{ram) 



1 

2^ 



2tt 



e i(m-n)9 de = g 



and 



oo 

') = ±- e in ^- e,) = 5(9 - 6') 



(33) 



(34) 



which can be easily derived from the closure relations 

c2tt 



^ />Z7T 

£ \n){n\=l, / d9\6)(9\ = l 



n— — oo 



and taking into account the definition 

1 / i ~ 



|n) 



'2it 



=e = (7i i 



(35) 



(36) 



We note in passing that both 9 and 6' are defined in 
the interval [0, 27r). In the coordinate representation, the 
e m6 jg an eigenfunction of L rcp = —ikd/d9 with corre- 
sponding eigenvalue nfi. This result can be expressed as 



(0|L(O)|n> = L Tep (e\n) = nh(9\n). (37) 
The wave function ijj now can be calculated: 

oo 

mt)= E (d\U(t)\n)(n\m) 

n— — oo 



/2vr 



If we define 



2?r 



(39) 



then the result of Eq. becomes that of Eq. l(T5|l . 



V. A LEWIS- RIESENFELD METHOD 

In this section, we briefly review the LR method and 
then apply it to solve the present problem. We shall show 
that LR method is not directly applicable, however, a 
simple modification concerning about the boundary con- 
dition makes it applicable to solving the problems with 
periodic boundary condition. 

Traditionally, to utilize the LR method j2^| solving a 
time-dependent system, we have to find an operator Q(t) 
such that 

ifi^=ih^ + {Q,H]=0,, (40) 
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and then find its eigenfunction ip\(9,t) satisfying 

Q(t)<p x (e,t) = \<p x (6,t), (41) 

with A being the corresponding eigenvalue. A wave func- 
tion ip\(0,t) satisfying Eq.QJ is then obtained via the 
relation 

fa(9,t) = e ia >®<px(9,t), (42) 
where a(t) is a function of time only, satisfying 

a x = <p x -\id/dt-H/h)< Px . (43) 
A general solution tp of Eq. (JJJ is then given by 

W,t) = J>(A)Va(M), (44) 



where g(X) is a weight function for A. 

To proceed let us assume the time-dependent invariant 
operator Q(t) takes the linear form 



Q{t) =a(t)L + b(t)6 + c(t). 



(45) 



in which a(t), b(t), and c(t) are time-dependent c-number 
functions to be determined. 

Substituting Eq. i|45|) into Eq. (|40() and solving these op- 
erator equations, we get 



but 



a(t) = a -—, b(t) = b , 



c{t) = c + b f n *^ H dT, 
Jo 1 



(46) 



(47) 



where ao, bo, and Co are arbitrary complex con- 
stants. Furthermore, substituting Eqs. H46J1 and l|47|) into 
Eq.ljlHJl, we find 

Q{t) = aoL{0) + bo9{0) + co = Q(0). (48) 

In other words, the invariant Q in the Heisenberg picture 
is precisely the linear combination of the initial canonical 
angular momentum L(0) and the initial azimuthal angle 
9(0) with an arbitrary constant cq. Note that in our 
system the L operator is also an invariant. 

How does the eigenvalue A evolve in time? Multiplying 
the factor e ta ^ on both sides of Ea. (|41|l . we get 



q(*)V>a(M = aVa(M). 



(49) 



Partially differentiating the both sides of Eq. l|49|l with 
respect to time and using Eq. (|40|l . we find 



A(t) = A(0), 



(50) 



To find a solution of Eq. , we have to solve Eq. i|41|) 
first. By solving Eq. l|41|l . we get 



y\(9,t) = exp 



where 



/*(*) 



A - c(t) 



K*) = -77T- 

a(t) 



Substituting Eq. (|S*T1) into Eq. lj4l?|) . we obtain 
'* [77 2 (t)+^(t)] 



a x (t) = a x (0) 



2Ih 



-dr. 



(51) 



(52) 



(53) 



where 



t](t) = h{t) - n 4 (r)/L (54) 
In deriving Ea. (|53|l we have used the two identities: 

/i= - , v = ~- (55) 

Here we see that in general a\(t) is a complex function. 

Although the form of ip\(9,t) = e iax(t) (f) X (9, t) is in- 
deed a solution of Eq.lQJ, however, it does not satisfy the 
periodic boundary condition [see Ea. (|12|) ]. This problem 
can be resolved by defining the total wave function ip(9, t) 
as the summation of all ipx {9 + 2nir, t) terms: 

oo 



Z exp 



n— — OQ 



iot\(t) + -fi(t)(6 + 2nic) 



2h 



v(t)(9 + 2nirf 



(56) 



It can also be transformed to the equivalent form below 
using the Poisson summation formula: 



y 2mv(t) 

OO 

z ex p 



exp 



ia x (t) + i^-9 2 c (t) + in 9 c (t) 



ih(n — no) 



2v(t) 



in (9-9 c (t)) 



(57) 



thus A is a constant. 



From these derivations it should be noted that when us- 
ing the LR method, the boundary conditions have to be 
carefully managed, otherwise one may get an incorrect 
result. 



VI. A COMPARISON TO VARIOUS 
APPROACHES 

In this section we shall show that Eq. Q56f) and 1571) 
can be cast into the forms of Eq.JTSJ and Eq. ij?T|) . To 
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proceed further, let us first borrow some parameters from 
Sec. II, and notice the simple result 

a(t)n h + b(t)9 c (t) + c(t) = a a n h + b 9 + c . (58) 

Comparing this result with Eq. (|5(J|I . we find that they 
are very similar. Now let us define 

X = a{t)n h + b(t)9 c {t) +c(t). (59) 

From Eq.JSHl), we get 

= n h + is{t)6 c {t), (60) 



rj(t) = l c (t) + v(t)9 c (t). 
Moreover, using the identity 

and the identity of v in Eq. (|55|l . we have 

, 2 = i? + 4K). 

Substituting Eq. (JHSJ) into Eq.(|S3I), we get 



e »A(0) exp , 



2Ti 



1 _ £ai 



(61) 



(62) 



(63) 



(64) 



Moreover, defining e laxty0 ^ and v$ as 



(65) 



we see clearly that Eqs. I|56f) and l|57|l become exactly the 
same as Eqs. (|18|) and (I21|l . These two approaches are 
thus verified to be equivalent. 

It is now interesting to discuss the physical meanings 
of l c (t) and 9 c (t) we have obtained. Although they are 
originated from the classical treatment, however, in what 
sense do they play a role of dynamic variables in the 
corresponding classical system should be further clarified. 
We would like to bring attention that in the Schrodinger 



picture within coordinate representation, the l c {t) is the 
expectation value of L — qRA(t) = —ihd/d9 — qRA(t). 
However, 9 c (t) is not the expectation value of 9 = 9 with 
respect to the wave function obtained in Eq. (|18|l , instead, 
it is merely the peak position of the wave packet [see 
Eq.|ZD]- 

In other words, the conventional Ehrcnfest Theorem 
is not directly applicable in this topologically nontrivial 
system. This consequence is due to the fact that we are 
not able to distinguish the phase between the angle 9 
and 9 + 2nir. Hence the 9 operator is not well-defined, 
only the e 10 is a well-defined operator, as has been 
demonstrated in Sec. IV. These facts cause the A los- 
ing its meaning as an expectation value of the Q operator. 

VII. SUMMARY 

In this article, we have studied the problem of a 
charged particle moving in a ring subject to a time- 
dependent flux threading it. After analyzing the prob- 
lem in a classical manner, various approaches includ- 
ing Fourier expansion method, time-evolution method, 
and Lewis-Riesenfeld method are considered and com- 
pared. In the Lewis-Riesenfeld approach, by appropri- 
ately managing the periodic boundary condition of the 
system, a time-dependent wave function can be obtained 
by using a non-Hermitian time-dependent linear invari- 
ant. The eigenfunction of the invariant can be realized as 
a Gaussian-type wave packet with the peak moving along 
the classical angular trajectory, while the distribution of 
the wave packet is determined by the ratio of the coeffi- 
cient of the initial angle to that of the initial canonical 
angular momentum. In this circular system, we find that 
although the classical trajectory and angular momentum 
can determine the motion of the wave packet; however, 
the peak position is no longer an expectation value of the 
angle operator, and the Ehrcnfest theorem is not directly 
applicable. 

Acknowledgment. — The authors are grateful to D. H. 
Lin and Y. M. Kao for discussion of the results. This 
work was supported partly by the National Science Coun- 
cil, the National Central University, and the National 
Center for Theoretical Sciences in Taiwan. 



[1] B. J. van Wees, Phys. Rev. Lett. 66, 2033 (1991). 

[2] C. Figueira de M. Faria, M. Dorr, and W. Sandner, Phys. 

Rev. A 55, 3961 (1996). 
[3] C. S. Tang and C. S. Chu, Phys. Rev. B 60, 1830, (1999). 
[4] M. T. Liu and C. S. Chu, Phys. Rev. B 61, 7645 (2000). 
[5] I. Guedes, Phys. Rev. A 63, 034102, (2001); H. Bekkar, 

F. Benamira, and M. Maamache, ibid. 68, 0160101 

(2003); I. Guedes, ibid. 68, 016102 (2003). 
[6] H. Maeda and T. F. Gallagher, Phys. Rev. Lett. 92, 

133004 (2004). 



[7] P. G. Luan and C. S. Tang, Phys. Rev. A 71, 014101 
(2005). 

[8] Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 
[9] A. G. Aronov and Yu. V. Sharvin, Rev. Mod. Phys. 59, 
755 (1987). 

[10] T. Chakraborty and P. Pietilainen, Phys. Rev. B 50, 8460 
(1994). 

[11] M. V. Berry, Proc. R. Soc. London A 392, 45 (1984). 
[12] Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 
(1987). 



6 



[13] R. Landauer and M. Buttiker, Phys. Rev. Lett. 54, 2049 
(1985). 

[14] Y. Gefen and D. J. Thouless, Phys. Rev. Lett. 59, 1752 
(1987). 

[15] K. Mullen, E. Ben.- Jacob, Y. Gefen, and Z. Schuss, Phys. 

Rev. Lett. 62, 2543 (1989). 
[16] M. Buttiker, Y. Imry, and R. Landauer, Phys. Lett. 96A, 

365 (1983); L. P. Levy, G. Dolan, J. Dunsmuir, and H. 

Bouchiat, Phys. Rev. Lett. 64, 2047 (1990). 
[17] V. Chandrasekhar et al., Phys. Rev. Lett. 67, 3578 

(1991); D. Mailly, C. Chapelier, and M. Benoit, Phys. 

Rev. Lett. 70, 2020 (1993); A. F. Morpurgo et al., Phys. 

Rev. Lett. 80, 1050 (1998). 



[18] R. Schuster et al., Nature (London) 385, 417 (1997); A. 

Fuhrer et al., Nature (London) 413, 822 (2001). 
[19] M. Bayer, M. Korkusinski, P. Hawrylak, T. Gutbrod, 

M. Michel, and A. Forchel, Phys. Rev. Lett. 90, 186801 

(2003). 

[20] A. Stern, Phys. Rev. Lett. 68, 1022 (1992). 
[21] J. Liu, B. Wu, and Q. Niu, Phys. Rev. Lett. 90, 170404 
(2003). 

[22] H. R. Lewis, Jr., Phys. Rev. Lett. 18, 510 (1967); J. 

Math. Phys. 9, 1967 (1968). 
[23] H. R. Lewis, Jr. and W. B. Riesenfeld, J. Math. Phys. 

10, 1458 (1969). 



